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1.  Introduction.  Quadratic  estimators  of  the  variance  components 
associated  with  random- effects  models  are  in  common  use.  If  the  experimental 
design  is  a  balanced  complete ly-nes ted  one,  if  the  usual  independence  and 
normality  assunptions  are  appropriate,  and  if  the  estimators  are  those 
derived  by  analysis -of- variance  techniques ,  then  each  of  the  estimators , 
save  that  of  the  residual  variance,  is  distributed  as  a  linear  difference 
between  two  independent  chi-square  variates.  More  generally,  any  quadratic 
form  in  multivariate-normal  data  and  in  particular  any  quadratic 
variance- component  estimator  based  on  such  data  is  either  distributed 
as  a  linear  difference  between  two  independent  chi-squares  or,  as  Press  [2] 
has  shown,  its  distribution  can  be  represented  as  a  mixture  of  such 
distributions. 

In  cases  where  a  quadratic  estimator  of  a  variance  component  can  take 
on  negative  values,  it  is  comnon  practice  to  replace  negative  estimates 
by  zero.  If  the  data  are  multivariate-normal,  the  distribution  of  the 
modified  estimator  is  identical  to  that  of  a  truncated  (at  zero)  linear 
difference  between  independent  chi-squares  or  is  a  mixture  of  such  truncated 
distributions. 

It  will  be  our  purpose  to  investigate  the  distribution  of  a  truncated 
linear  difference  between  independent  chi-square  variates.  Rather  than 
assume  that  the  truncation  is  from  below  and  at  zero,  we  will  consider  the 
more  general  situation  where  the  truncation  is  either  from  above  or  below 
or  both  and  at  arbitrary  'values.  In  particular,  we  will  obtain  expressions 
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for  the  probability  density  function,  distribution  function,  moments, 
moment  generating  function,  and  characteristic  function  of  such  a  distribution. 
In  the  case  of  the  probability  density  function,  distribution  function, 
and  moments,  emphasis  will  be  on  obtaining  representations  that  appear  to 
be  convenient  computational  forms.  Extensive  use  will  be  made  of  known 
results  on  hypergeometric  functions. 

Previously,  Press  [2]  obtained  the  probability  density  function  of 
a  linear  difference  between  independent  chi-squares  in  terms  of  a  confluent 
hypergeometric  function  of  the  second  kind.  Wang  [5]  found  a  simpler 
expression  for  the  probability  density  function  for  the  special  case 
where  the  degrees  of  freedom  associated  with  both  chi-squares  are  even. 

She  also  gave,  for  that  same  special  case,  convenient  expressions  for 
the  moments  of  that  distribution,  after  its  truncation  from  below  at  zero. 

2.  Notation  and  preliminary  results^  Define 

X  -  aWx  -  BW2, 

where  a  and  B  are  positive  constants,  and  and  W2  are  indpendent  chi-square 
variates  having  m  and  n  degrees  of  freedom,  respectively.  Truncation  of 
X  from  above  and/or  below  yields  the  random  variable 
Y  ■  yj  if  X  >  y^f 

*  X  ify1>X>y0, 

“  y0  if  x  <  y0> 

where  yQ  and  y^  are  known  constants  with  -»  <  yQ  <  y^  <  ®. 

If 

A  chi-square  variate  with  r  degrees  of  freedom  has  kth  moment  2  (r/2)^, 
where  (a)Q  -  1  and  (a)k  «  a(a+l) . . . (a+k-1) ,  k  »  1,  2,  ...;  moment 
generating  function 

M(0)  -  (1  -  20)'r/2, 
and  characteristic  function 

C(0)  -  (1  -  2i0)'r/2. 


0  <  1/2; 
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It  follows  that  X  has  kth  moment 


U)  (2«)kIj.,  (•)  (-03  (n/2)j  (m/2)k.r 

where  c  *  8/a;  moment  generating  function 

(2)  M^(8)  -  (1  -  2a9)'m/2  (1  ♦  289)'n/2,  (-28)' 


(~28) ~ 1  <  e  <  (2a)'1; 


and  characteristic  function 


C^(0)  -  (1  -  2ai6)“m/2  (1  ♦  28i6)'n/2. 


U(c,  d;  x)  -  (l/r(c)]  fl  e‘xt  tc‘X  (1  ♦  t)d'c_1  dt,  c  >  0,  x  >  0, 
to  be  a  confluent  hypergeometric  function  of  the  second  kind,  and  put 

m  Hr  -  .  h-  tl  -  t(r*s'2)/Z  e't/(2a)  Ufs/2,  (r*s)/2;  (a»b)t/(2ab)l 

wJ  g(r,  s,  a,  o,  t)  * - ?rJiT75 — r72 — <75 - 

2ir*sj/^  ^r/i  bs/z  r(r/2) 

for  r  -  1,  2,  ... ;  s  •  1,  2,  ... ;  a  >  0,  b  >  0;  and  t  >  0. 

From  (2. ,  we  have  that  the  probability  density  function  of  X  is 


fx00  ■  g(m,  n,  a,  8;  x) 

»  g(n,  m,  8,  a;  -x) 

.  ^ 2  r(m/2  ♦  n/2  -  1) _ 

2(m+n)/2  g  [(?+i)/2] (m+n-^)/2  r(m/2)  T(n/2) 


if  x  >  0, 
if  x  <  0, 


if  x  -  0. 


It  will  prove  convenient  to  express  the  distribution  functions  of 
X  and  Y  and  the  moments,  moment  generating  function,  and  characteristic 
function  of  Y  in  terms  of  the  function 

G(r,  s,  a,  b;  k,  z,  u)  -  /“  tk  ezt  g(r,  s,  a,  b;  t)  dt  for  u  >  0, 


for  u  »  0, 


where  k  ■  0,  1,  2,  ... 


Note  that,  for  v  <  0, 

(6)  /y  tk  ezt  g(s,  r,  b,  a;  -t)  dt  -  (-l)k  G(s,  r,  b,  a;  k,  -z,  -v). 

The  distribution  function  of  X  is 

Fx(x)  -  G(m,  n,  a,  8;  0,  0,  x)  +  G(n,  m,  8,  a;  0,  0,  »)  if  x  >  0, 

«  G(n,  m,  8,  a;  0,  0,  ®)  -  G(n,  m,  8,  a;  0,  0,  -x)  if  x  <  0. 
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Denote  by  P_  («)  the  distribution  function  of  an  F  random  variable 
r 

having  r  and  s  degrees  of  freedom.  Note  that 

Pm,n^ (nS)/ Cmcx) }  -  Fx(0)  -  G(n,  m,  8,  o;  0,  0,  «). 

The  distribution  function  of  Y  is 
Fy(y)  -  i  ify>yx, 

-  Fx(y)  if  y0  <  y  <  yv 

-  0  if  y  <  yQ. 

If  yQ  and  y^  are  finite,  Fy(*)  is  discontinuous  at  those  points.  It 
makes  a  junp  of  Fx(yQ)  at  yQ  and  1  -  Fx(y^)  at  y^.  At  other  points,  the 
probability  density  function  of  Y  is 

Fy(y)  -  fx(y)  ify0<y<y1, 

•  0  if  y  <  yQ  or  y  >  yj. 

Define 

H(ur  tu:  k,  zl  «  Gfm.  n.  a,  8;  k,  z,  Uj)  -  G(m,  n,  a,  8;  k,  z,  Uq) 

if  0  <  Uq  <  Uj, 

G(«»  n»  a»  6;  «x)  ♦  (-Dk  G(n,  m,  8,  a;  k,  -z,  -Uq) 

if  Uq  <  0  <  ur 

■  (-1)  [G(n,  m,  6,  a;  k,  -z,  -Uq)  -  G(n,  m,  8,  o;  k,  -z,  -Uj)] 

if  uQ  <  ux  <  0; 
and 

H*(up  Uq;  k,  z)  -  H(ur  UqJ  k,  z)  ♦  u*  eZU9  Fx(uQ)  ♦  u\  ezui  [1  -  F^)] 

for  -•  <  Uq  <  ^  <  •», 

Uq;  k,  z)  -  H(»,  Uq;  k,  z)  ♦  uj  ezu«  Fx(uQ)  for  —  <  Uq, 

H*(Uj ,  k,  z)  »  H(Up  k,  z)  ♦  uj  ezui  [1  -  Fx(u1)]  for  ^  <  •, 
and 

H*(»,  k,  z)  ■  H(«,  k,  z) , 
with  0°  *  1. 

The  kth  moment  of  Y  is  H*(yp  ygj  k,  0),  and  its  moment  generating 
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function  and  characteristic  function  are 

My(6)  -  H*(yr  y0;  0,  0),  0  >  (-20)' 1  if  yQ  <  0  and  0  <  (2a)'1  if  y2  >  0, 

and 

Cy(6)  -  H*(yp  yQ;  0,  i0), 
respectively. 

It  is  clear  that  the  problem  of  evaluating  the  probability  density 
function  of  X  or  Y  can  be  reduced  to  one  of  evaluating  g(r,  s,  a,  b;  t) 
for  appropriate  values  of  that  function’s  parameters.  Moreover,  to 
evaluate  the  distribution  function  of  X  or  Y  and  the  moments  and  moment 
generating  function  of  Y,  it  suffices  to  evaluate  G(r,  s,  a,  b;  k,  6,  u) 
for  appropriate  values  of  r,  s,  a,  and  b  and  for  k  ■  0,  1,  2,  . ..,  9  real, 
and  u  >  0.  For  the  most  part,  we  will  find  it  convenient  to  carry  out 
the  analysis  for  all  such  (k,  9,  u)  triplets,  rather  than  to  separately 
consider  only  those  particular  triplets  that  are  relevant  to  the  evaluation 
of  the  distribution  function,  the  moments,  and  the  moment  generating  function, 
respectively.  Because,  in  (5),  t  e  is  'absorbed'  by  g(r,  s,  a,  b;  t), 
the  more  general  approach  is  no  more  difficult.  It  has  the  advantage  of 
avoiding  repetition.  Since  the  expressions  for  both  the  moment  generating 
function  and  the  moments  are  to  be  obtained  directly,  the  derivatives  at 
zero  of  expressions  for  the  moment  generating  function  can  be  compared 
with  the  corresponding  expressions  for  the  moments ,  so  as  to  obtain  a 
check  on  the  correctness  of  the  results. 

Expressions  for  the  characteristic  function  Cy(e)  of  Y  will  not  be 
given  explicitly;  however,  substitution  of  ie  for  e,  in  any  of  the 
expressions  indicated  for  the  moment  generating  function  My  (8) ,  produces 
a  valid  representation  for  Cy(6),  provided  that  representation  is 
meaningful . 


3.^  Evaluation  of  g(r,  s,  a,  b;  x).  Us^  ol  relationships.  It  is 
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easy  to  show  that,  for  c  >  0, 

(7)  g(r,  s,  a,  b;  t)  ■  c  g(r,  s,  ca,  cb;  ct). 

The  recurrence  relations  for  confluent  hypergeometric  functions  given 
by  (13.4.15)  -  (13.4.20)  of  Slater  (4)  yield 

(8)  g(r,  s,  a,  b;  t)  *  (b/a)(r-2)'1  {(t(aM+b"1)  ♦  s  -  r  ♦  4]  g(r-2,  s*2,  a,  b;  t) 

♦  (b/a) (s+2)  g(r-4,  s+4,  a,  b;  t)}, 

(9)  g(r,  s,  a,  b;  t)  ■  b(a^b)'1 (r-2)'1  {{t(a‘l+b'1)  ♦  s  ♦  r  -  4]  g(r-2,  s,  a,  b;  t) 

-  (t/a)  g(r-4,  s,  a,  b;  t)}, 

(10)  g(r,  s,  a,  b;  t)  -  [a(r-2)]M  (bs  g(r-2,  s+2,  a,  b;  t)  ♦  t  g(r-2,  s,  a,  b;  t)J, 

(11)  g(r,  s,  a,  b;  t)  -  (a+b)'1  (b  g(r-2,  s,  a,  b;  t)  ♦  a  g(r,  s-2,  a,  b;  t)J, 

(12)  g(r,  s,  a,  b;  t)  «  b(a*b)'l(r-2)‘l  {(t(a'l+b‘l)  ♦  s]  g(r-2,  s,  a,  b;  t) 

♦  (bs/a)  g(r-4,  s+2,  a,  b;  t)}, 


and 

(13)  g(r,  s,  a,  b;  t) 


(b/a) (r-2)”1 


{[t(a"l*b”1)  *  s ] 


fT  (  **-  ? 

e>  v*  -  t 


c**)  i  K* 
w  -»  *■'» 


-  (t/a)  g(r-4,  s+2,  a,  b;  t)}. 

Hie  probability  density  function  of  a  chi-square  variate  with  r 
degrees  of  freedom  is 

Pr(t)  -  [2 r/2  T(r/2) ]’ 1  t^2*'1  e't/2,  0  <  t  <  -, 

•  0,  elsewhere. 

The  domain  of  the  function  g  can  be  extended  by  putting 
g(r,  0,  a,  b;  t)  =  (1/a)  pr(t/a) 

for  r  *  1,  2,  ...;  a>0;  b>0;  t>0.  The  relationships  (7)  -  (13) 

remain  valid.  Thus,  tables  of  the  chi-square  density  may  be  of  use  in 
evaluating  the  densities  of  X  and  Y. 

Case  (i):  r  an  even  integer.  For  c  >  0  and  p  =  1 ,  2 ,  . . . , 

U(c,  c+p;  x)  ■  [l/r(c)l  jr  j  (P:1)  /“  e"xt  dt 

(14)  -  if!  pj1)  (c)j  x'c'i. 


Substitution  of  this  expression  into  (3)  gives 
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(15)  g(r,  s,  a.  b;  x)  -  [a  2r/2  r(r/2))'‘  [a/(a*b))s/2  e'x/(2a) 

•  (T/f *)  (s/2)j  (2b/(a.b)]j  (x/a)r/2'j_1, 

for  r  *  2,  4,  6,  ... 

Case  (ii):  s  an  even  integer.  Denote  by  Pu(*)  the  distribution 
function  of  a  chi-square  variate  with  u  degrees  of  freedom.  For  p  »  1,  2,  ... 
and  2c  8  1,  2,  ...» 

U(p,  p+c;  x)  -  [l/r(p)J  rx  e”x^v  l)  (v  -  l)?'1  v -c_1  dv 

-  U/TCp)]  eX  l]: ]  P:1)  C-DP^'1  K  e_XV  vC+^2  dv 

-  [i/r(p)l  ex  p:1)  (-i)P*^1  x-°-j  r(c*j) 

*  *  P2c+2j^2x^* 

Substituting  this  expression  into  (3),  we  find 

g(r,  s,  a,  b;  x)  -  [b  2s/2  r(s/2)f 1  (b/(a*b)]r/2  ex/(2b) 


•  I*£?' 1  f57-'1)  Cr/2)^  [2a/(a»b))j  (-x/b)s/2'j_1 

'  !l  '  "r.2j!x(a*b)/tab)”- 

for  s  «  2,  4,  6,  ... 

Case  (iii):  Both  r  and  s  are  odd  integers.  Let 


(u). 


cj(u-  v)  ’  TO* 


j  »  1,  2,  ... 


Take 


*(u,  v;  z)  «  Cj(u,  v)  z) 

to  be  a  confluent  hypergeometric  function  of  the  first  kind.  This  series 
is  absolutely  convergent  for  all  values  of  u,  v,  and  z,  save  v  *  0,  -1,  -2,  ... 
Its  remainder  after  q  terms  will  be  denoted  by  the  symbol 
Rq(u,  v;  z)  =  I",q  Cj  (u,  v)  zK 
The  numerical  evaluation  of  U[s/2,  (r+s)/2;  (a+b)x/(2ab) J  and 
consequently  of  g(r,  s,  a,  b;  x)  can  be  accomplished  for  r  *  1,  3,  5, 
s  *  1,  3,  5,  ...,  by  applying  the  well-known  relationship 
(16)  U(v,  1+p;  t)  «  ((-1)P  p!  r(v-p)JM  L(v,  1+p;  t) 

‘  [(p-i)!/r(v))  Cj(v-p,  1-p) 
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for  p  ■  0,  1,  2,  ...;  where,  taking  ^(t)  *  r»(t)/r(t)  to  be  the  logarithmic 
derivative  of  the  gamma  function  and  y  *  -  ^(1)  3  0.57721  56649...  to  be 
Euler's  constant,  L(v,  1+p;  t)  is  obtained  by  adding  together  the  three 


convergent  series 


(17)  (*Y  -  In  0  *(v,  1+p;  t), 

(18)  c^(v,  1+p)  tJ  [^(1+j+p)  *  <Hv+j)], 
and 


(19)  IJ«0  cj(y*  1+P)  tj  lL  u*l‘ 

Here,  the  meaningless  sums  L*  and  j[®  are  interpreted  as  zero. 

J  U 

The  series  (17) ,  (18) ,  and  (19)  can  be  approximated  by  their  partial 


sums.  We  now'  show,  for  0  <  v  <  1  ♦  p  and  t  >  0,  how  to  estimate  the 


errors  incurred  in  these  approximations.  These  results,  together  with 
(16)  and  tables  of  the  gamma  function,  allow  us  to  calculate  U(v,  1+p;  t), 
for  such  v.  p,  and  t  values,  with  known  accuracv. 

Denote  by  [v]  the  largest  integer  that  is  less  than  v.  By  assumption 
0  <  v  <  1  +  p  and  t  >  0,  so  that 

(20)  0  <  Rq(v,  1+p;  t)  <  R.q(l+(v],  1+p;  t). 

Furthermore,  since  u'(y)  >  0  and  v' ' (y)  <  0  for  y  >  0,  we  have 

0  <  Cj (v,  1+p)  t^  [<J/(l+j+p)  -  ’Kv+j)] 

(21)  <  Wl+q+p)  *  'Hv+q)]  Rq(v,  1+p;  t) 

(22)  <  fi|>(l+q+p)  -  Kv+q)]  Rq(l+[v],  1+p;  t). 

Also,  for  q  >  1, 

0  <  £.q 

<  ["=q  j  Cj (V,  1+p)  ^ 

3  vt  (1+p)'1  Rq.jfl+v,  2+p;  t) 

(23)  <  vt  (1+p)'1  Rq.i(2+[v],  2+p;  t). 


The  recurrence  relationship 
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(24)  w(w+l)  4>(y,  w;  t)  -  (w+l)(w+t)  4(y,  w+1;  t) 

♦  t(w-y+l)  4(y,  w+2;  t)  =  0; 

together  with  the  formulas 

*(j,  j;  t)  -  e1 
and 

*(j»  l+j;  t)  *  (-1)^  ji  [1  -  el  Cj.jC-t)],  j  ■  1,  2,  ...» 
where  e  (t)  «  K.#  tJ/j  ! ;  can  be  used  to  construct  convenient  procedures 

»  J 

for  computing  R  (l+[v] ,  1+p;  t)  and  R  ,(2+[v],  2+p;  t)  and  thus,  in  light 
4  4  A 

of  (20),  (22),  and  (23),  can  be  used  for  estimating  errors  resulting  from 
the  use  of  partial  sums  in  place  of  the  series  (17) ,  (18) ,  and  (19) ,  and 
for  determining  the  number  of  terms  that  need  be  included  in  each  partial 
sun. 

For  many  v,  p,  and  t  values,  it  may  be  advantageous  to  modify  the 
outlined  procedures  for  approximating  the  series  (17) ,  (18) ,  and  (19)  with 
known  accuracy.  Existing  tables,  Kurmer's  transformations ,  recurrence 
relationships  like  (24),  and/or  other  techniques  fran  the  theory  of 
confluent  hypergeometric  functions  (see  e.g.  [3]  or  [4])  can  often  be 
utilized  to  approximate  $(v,  1+p;  t)  and  consequently  the  series  (17) 
in  more  efficient  fashion  than  the  partial-sum  approach  described  above. 

Our  more  efficient  approximation  for  4>(v,  1+p;  t)  also  yields  R  (v,  1+p;  t) 
which,  upon  applying  (21),  can  also  be  used  to  estimate  the  error 
resulting  from  approximating  (18)  by  one  of  its  partial  sums.  Note  that 
this  error  estimate  is  more  precise  than  the  one  based  on  R^(l+[v],  1+p;  t). 
Similarly,  it  may  be  possible  to  approximate  <t(l+v,  2+p;  t)  and  consequently 
Rq.jU+v,  2+p;  t)  to  known  accuracy  with  techniques  not  requiring  the 
computation  of  R  ,(2+[v],  2+p;  t) ;  thus  permitting  more  precise,  and  possibly 

4  A 

easier,  estimation  of  errors  resulting  from  truncation  of  the  series  (19). 
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4.^  Evaluation  of  G(r,s,  a,  b;k,e,x).  Useful  relationships. 

Using  (7),  we  find 

G(r,  s,  a,  b;  k,  0,  x)  «  c  G(r,  s,  ca,  cb;  k,  8/c,  cx). 

Upon  substituting  the  expression  (8)  for  g  into  the  definition  (5), 
we  obtain 

G(r,  s,  a,  b;  k,  6,  x) 

-  (b/a) (r-2)'1  {(a'^b'1)  G(r-2,  s+2,  a,  b;  k+1,  9,  x) 

♦  (s-r*4)  G(r-2,  s+2,  a,  b;  k,  9,  x) 

♦  (b/a)(s+2)  G(r-4,  s+4,  a,  b;  k,  8,  x)  . 
Recurrence  relations  for  G,  corresponding  to  (9)  -  (13),  can  easily  be 
obtained  in  the  same  way. 

Hie  domain  of  G  can  be  extended  by  defining  G(r,  0,  a,  b;  k,  0,  x) 
in  terns  of  g(r,  0,  a,  b;  •)  through  application  of  (5).  Our  recurrence 
relations  for  G  remain  valid.  Thus,  tables  of  the  chi-sauare  distribution 
function  may  prove  useful  in  evaluating  the  distribution  functions  of  X 
and  Y  and  the  moments  of  Y. 

By  making  use  of  (1),  (2),  (4),  and  (6),  we  obtain 

(25)  G(r,  s,  a,  b;  k,  0,  «)  •  {-l)k+1  G(s,  r,  b,  a;  k,  0,  ”) 

♦  (2a)k  (k)  (-b/a)j  (s/2).  (r/2)fc_j 

and 

(26)  G(r,  s,  a,  b;  0,  0,  •)«(!-  2a0)’r/2  (1  +  2b9)*s/2 

-  G(s,  r,  b,  a;  0,  -8,  »), 

-l/(2b)  <  8  <  l/(2a). 

Case  (i):  r  an  even  integer.  Using  (15),  we  find 
G(r,  s,  a,  b;  k,  0,  x) 

-  [a2r/2  r(r/2)J*‘  [a/(a+b)]s/2 

*  l]i]~ 1  (s/2)j  [2b/ (a+b) p  a*r/2+:i+1 

*  !*  tr/2+k'j'!  exp[-t(l-2a0)/ (2a) j  dt 


11 


-  [2r/2  T(r/2)]'1  [a/(a*b)]s/2  ak 

*  r/Jt'1  (r/fX)  (s/2)j  (2b/(a*b))^  (l-2a8)'r/2'k+j 
.  |X(l-2a9)/a  yr/2+k-j-l  #-y/2  ^ 

k  s/2  pr/2-1  (r/2-j)k  ^  .  j 

(27)  -  (2a)*  [a/(a*b))s/Z 

•  (l-2a8)'r^Z’k^  Pr*2k-2)[x(l‘2*e)/*,; 
for  r  ■  2,  4,  6,  6  <  l/(2a);  0  <  x  <  *. 


<  In  particular, 

(2b)  G(r,  s,  a,  b;  k,  8,  -)  -  (2a)k  [a/(a*b))s/2 

•  &)’  d-2ae)'r/2'k*3: 

for  r  ■  2,  4,  6,  0  <  l/(2a). 

Applying  the  well-known  relationship 

P2v(2x)  -  [r(vl)]“  xv  4(v,  1+v;  -x) 
to  (27),  we  obtain,  for  0  <  x  <  •,  the  alternate  representation 
G(r,  s,  a,  b;  k,  0,  x) 

-  xk  [x/(2a)]r/2  [a/(a*b)]s/2 

.  vr/2-1  (5/2)j  r  2ab  J 

2j-o  3T  IrTZ-j-DV  (r/MU  IxfiFE? 

•  ♦[r/2*k-j,  r/2+k-j+l;  -x(l-2a0)/(2a)J ; 

r  •  2,  4,  6,  8  <  l/(2a). 

Case  (ii):  s  an  even  integer.  By  working  with  the  joint  probability 
distribution  of  the  random  variables  Wj  and  X,  we  obtain  the  representation 
(29)  fx(t)  *  (1/8)  r^,t/a)  p„(u)  Pn[(ou-t)/6)  du 

for  the  probability  density  function  of  X. 

Using  (29),  together  with  (4)  and  (6),  we  find;  for  s  ■  2,  4,  6, 
-l/(2b)  <  8  <  l/(2a);  0  <  x  <  ®; 


G(r,  s,  a,  b;  k,  0,  x) 

-  (-l)k+1  G(s,  r,  b,  a;  k,  -0,  ®) 
♦  (1/b)  /*  dt  tk  e0t  T 


max(0,t/a) 


p^(u)  ps((au-t)/b]  du 


12 


-  i-l)k+1  G(s,  r,  b,  a;  k,  -6,  «) 

♦  ak  (a/b)  !xi3  dv  vk  ea6v  f^xio.v)  Prfu)  Psla(U'v)/b]  du 

-  (-I)11*1  G(s,  r,  b,  a;  k,  -8,  -) 

♦  Gj(r,  s,  a,  b;  k,  8,  x)  ♦  G2(r,  s,  a,  b;  k,  0,  x) 


where 


G^r,  s,  a,  b;  k,  6,  x) 

-  ak  (a/b)  f*/a  du  pr(u)  /?„  vk  ea°v  ps[a(u-v)/b]  dv 

»  ak  f*/a  du  pr(u)  K  (u  -  by/a)k  ea0tu-by/a)  ps(y)  dy 
.  ak  r/a  du  pr(u)  ea6u  Ik.t  (k)  u>  (-b/a)k'j  £  P,M  dy 

-  (-2b)k  (l-2a0)"r/2  (l+2b8)"s/2’k 

•  l)  .  (■)  (r/2L  Cs/2)v_.  (-a/b)3  t(l*2b8)/(l-2a8)]j  P  2.[x(X-2ae)/a) 

J  j  *  j  ^ 

and 

G2(r,  s,  a,  b;  k,  6,  x) 

*  ak  (a/b)  ,^/a  du  pr(u)  !**  vk  ea6v  p$[a(u-v)/b]  dv 

-  ak  (a/b)^2  [2s'2  rcs/Hl"  du  pr(u)  e'au/'2b> 

.  p/2-1  (S/2-1)  us/2-j-l  c.x)j  jx/a  yk*j  eav(l*2b8)/(2b)  dy 

-  (-2b)^  (b/(a*b)lr/z  (l»2b8)'5/2'k  el(1'’2b6)/(2b) 

•  ip1  (^}ja*2be)j 

es/2*k-i-l!-x(1*2be)/(2b)l  {1  ‘  Pr*2j  f*(a»b)/(ab) I >. 

Note  that  G(s,  r,  b,  a;  k,  -0,  »),  which  appears  on  the  right  hand 

side  of  (30),  can  be  evaluated  from  (28).  Expression  (28),  together  with 
(25)  and  (26),  can  also  be  used  to  evaluate  G(r,  s,  a,  b;  k,  0,  »)  and 
G(r,  s,  a,  b;  0,  9,  «)  for  s  -  2,  4,  6,  ... 


Case  (iii):  General  case.  Take 


F(a,  b,  c;  z) 


J1  tb'1  (l-t)c'b_1  (l-tz)'a  dt, 

*  0 

Re  c  >  Re  b  >  0,  |arg(l-z)|  <  ir, 


to  be  the  analytic  continuation  of  the  hypergeometric  function. 
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For  the  Laplace  transform  of  yv' 1  U(c,  d;  y),  we  have 
ft  e"2^  yv  l  U(c,  d;  y)  dy 

"  ^'T{lVciv'‘d)d^  ^  F(v’  1+v'd>  1+c+v'd;  1_z) 

providing  Re  v  >  0,  Re  d  <  Re  v  ♦  1,  Re  z  >  0  (equation  (3. 2. SI)  of 
Slater  [3]  and  its  extension  by  analytic  continuation);  so  that 
G(r,  s,  a,  b;  k,  0,  “) 

■  [2(m)/2  ar/2  bs/2  r(r/2)f 1  (2ab/(a+b)]r/2+s/2+k 

*  ro  yr^2*k'1  exp(-yb(l"2a6)/ (a+b)]  U(s/2,  r/2+s/2;  y)  dy 

.  f  a  is/ 2  (  b  ) r/2+k  T(k+1)  r(r/2+s/2+k) 

Uaj  1*5*  T(r7'2T’T  (s/2+k+ 1) 

•  F[r/2+s/2*k,  k+1,  s/2+k+l;  a(l+2b0)/(a+b)] 

for  0  <  l/(2a). 

By  making  use  of  linear  transformation  formulas  for  the  hypergeometric 
function  ((9.5.2)  and  (9.5.3)  in  Lebedev's  book  [1]),  we  obtain  alternate 
representations 

G(rt  s,  a,  b;  k,  0,  -) 

-  tM*  r_ajs/2  (  b  ^r/2-1  r(k*l)  r(r/2+s/2+k)  „  ,_.„-k-l 
(x n  UaJ  W  OT  T(r77)’  :(s/7+k*l)  (l'4a0) 


•  F[l-r/2,  k+1,  s/2+k+l;  - (a/b) (l+2b9)/(l-2a0) j , 


0  <  l/(2a), 


G(r,  s,  a,  b;  k,  0,  «) 

,  O.^k  r  a  is/2  r(k+l)  r(r/2+s/2+k)  „  ._„,-r/2-k 
lz*J  w  r (r/2)  r(s/2»knr  (1'2a9J 

•  F[l-r/2,  s/2,  s/2+k+i;  a(l+2b6)/(a+b) J , 

In  particular,  the  representation  (31)  yields 
G(r,  s,  a,  b;  k,  0,  «°) 

.  f2a)k  (-*  )s/2  fA.)r/2-l  r(k+l)  r(r/2>s/M) 

1  J  W*  r (r/2)“ r  (s/2+k+lj 

*  F(l-r/2,  k+1,  s/2+k+l;  -a/b). 


0  <  1/ (2a) . 


(33) 
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This  expression  is  an  especially  fruitful  computational  form. 

For  r  *  2,  4,  6,  . ..,  it  reduces  to 

(34)  G(r,  s,  a,  b;  k,  0,  •)  ■  (2a)k  (^fg) r/2+s/2‘1 

rr/2-i  (1-r/2-s/2-k).  (r/2-j),  . 

*  ifJ  1  - JT2 - -  (W. 

providing  an  alternative  to  (28)  for  computing  purposes. 

For  a/b  <  1  and  r  »  1,  3,  5,  ...,  terms  (r-l)/2,  (r+l)/2,  (r+3)/2,  ... 
of  the  right  hand  side  of 

(l-r/2) .  (k+1).  - 

(35)  F(l-r/2,  k+1,  s/2+k+l;  -a/b)  -  £je0  (*a/b)J 

comprise  a  convergent  alternating  series,  whose  remainder  is  smaller  in 
absolute  value  than  the  first  neglected  term  and  has  the  same  sign;  so 
that,  by  making  use  of  the  series  representation  (35),  G(r,  s,  a,  b;  k,  0,  •) 
can  be  calculated  to  any  desired  accuracy. 

In  evaluating  G(r,  s,  a,  b;  k,  0,  °°)  from  (25)  for  s  *  2,  4,  6,  . . . , 
formula  (34)  can  be  used  as  cn  alternative  to  (28)  for  computing 
G(s,  r,  b,  a;  k,  0,  “).  Similarly,  if  a/b  >  1  and  s  *  1,  3,  5,  ..., 

G(s,  r,  b,  a;  k,  0,  °°),  which  appears  on  the  right  hand  side  of  the 
expression  (25)  for  G(r,  s,  a,  b;  k,  0,  ®),  can  be  evaluated  by  making 
use  of  (31)  and  the  series  representation  (35). 

From  (32),  we  obtain 

G(r,  s,  a,  b;  k,  0,  «)  *  (2a)k  (r/2+s/2)k  £k-0  (k)  ((s/2) ./(r/2+s/2) .] 

•  (-lp  (l+b/ap  (1  -  Pr>s+2j((b/a)(s+2j)/r]}, 
so  that  this  function  can  also  be  computed  from  tables  of  the  F  distribution 
function  or  the  incomplete  beta  function. 

For  0  <  x  <  »,  we  find 

2(r+s)/2  ar/2  bs/2  r(r/2)  r(s/2)  G(r>  s>  a>  b;  k>  0f  x) 

-  ft  du  us/2_1  (l*u)r/2-1  /*  tr/2+s/2+k_1  exp{-t[u(l+a/b)+l-2a8]/(2a) }  dt 
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-  (r/2+s/2+k) ’ 1  xr/2+s/2+k  exp[-x(l-2a0)/(2a)] 

•  (lm)r^21  exp[-ux(a+b)/(2ab)} 

•  *{1,  r/2+s/2+k+l;  x[u(l+a/b)+l-2a0]/(2a)}  du 

-  (r/2+s/2+k) * 1  xr/2+s/2+k  cxp[-x(l-2a0)/(2a)] 


•  T(s/2+j)  U[s/2+j,  r/2+s/2+j ;  x(a+b)/(2ab)] 


-  xr/2+s/2+k  exp[-x(l-2a0)/(2a)] 


(36)  *  (r/zis/2^}—  (x(a+b)/(2ab)]^  U[s/2*j,  r/2*s/2*j;  x(a*b)/(2ab)] 

•  ♦[j+l,  r/2+s/2+k+j+l;  x(l-2a0)/(2a)] , 

0  <  l/(2a). 


«  x 


Using  Ki -Timer's  transformation,  we  obtain  an  alternate  representation 
2(r*s)/2  ar/2  bs/2  r(r/2)  r(s/2)  G(r,  s,  a,  b;  k,  0,  x) 

'r/2"/2*k 

•  U[s/2*j,  r/2+s/2+j;  x(a+b)/(2ab)J 

•  *[r/2+s/2+k,  r/2+s/2+k+j+l;  -x(l-2a0)/(2a)j , 

0  <  x  <  08 ,  0  <  1/ (2a) . 


We  now  consider;  for  r  *  1,  3,  5,  s  *  1,  3,  5,  estimation 

of  the  error  incurred  in  approximating  the  series  appearing  on  the  right 
hand  side  of  (36)  by  a  partial  sum. 

We  will  need  the  simple  inequalities 

0  <  r(c)  U(c,  d;  t)  <  r(c-6)  U(c-6,  d;  t), 

t  >  0,  c  >  c-6  >  0, 
and 

0  <  *(c,  d;  t)  <  e*, 


t  >  0,  d  >  c  >  0, 


and  formula  (14). 
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Putting  r*  ■  r  ♦  1  and  s*  ■  s  -  1,  we  find  that,  for  r  -  1,  3,  5, 
s  ■  1,  3,  S,  and  for  any  positive  integer  v, 

0  <  [r(s/2*j)/(r/2+s/2+k)j+j] [x(a+b)/(2ab)p 

•  U[s/2«-j ,  r/2+s/2+j;  x(a+b)/(2ab)]  $>[j-*-l,  r/2+s/2+k+j+l ;  x(l-2ae)/(2a) ) 

(37)  <  £v  lr(sV2+j)/(rV2^sV2+k)j+1][x(a^b)/(2ab))j 

•  U(s*/2*j,  r*/2+s*/2*j;  x(a+b)/(2ab)J 

*  rV2+s*/2+k+j+l;  x(l-2a6)/(2a)J 

<  ex(l-2a0)/(2a)  ^  [r(s*/2+j)/(r*/2*sV2+k)  .+1J [x(a*b)/(2ab)] * 

•  U[s*/2*j ,  r*/2+s*/2+j ;  x(a+b)/(2ab) j 
•  [x(a+b)/(2ab))’s^2  ex(1_2a0)/(2a) 

‘  fjuy  f  (rV2+s*/2+k)  j+1  J  * 1  [p^2*1  f'*1)  r(s*/2+p+j) 

[x(a+b)/(2ab)fp 

(38)  "  tx(a+b)/(2ab)fs*/2  ex(1‘2a6)/(2a) 

•  ^p2'1)  (x(a*b)/(2ab)fp 

*  ^Tr*/2+s */ 2+k) f r V 2+k-p)  *  %l\ 

0  <  x  <  ■>,  0  <  l/(2a) . 

For  r  ■  1,  3,  S,  s  *  3,  5,  7,  we  note  that  the  error 
bound  (37)  equals 

2(r+s)/2  ar*/2  fas*/2  r(r*/2)  r(s*/2)  x“r/2‘s/2*k 

(39)  ex(l-2a0)/(2a)  G(r*.  s*.  a,  b;  k,  0,  x) 

'  £j«J  [r(s*/2+j)/(r/2+s/2+k) j+1J [x(a+b)/(2ab) 

•  U(sV2+j,  r/2+s/2+j ;  x(a*b)/(2ab)] 

•  ♦[j+1,  r/2+s/2+k+j*l;  x(l-2a0)/(2a)J. 

Here,  G(r*,  s*,  a,  b;  k,  0,  x)  can  be  computed,  for  example,  from  (27); 

U[s*/2+j ,  r/2+s/2+j;  x(a+b)/(2ab) j  can  be  calculated  from  (14);  and 
♦[)'*'!»  r/2+s/2+k+j+l;  x(l-2a9)/(2a)]  can  be  evaluated  by  using  the  recurrence 
relationship  (24)  together  with  its  two  accompanying  formulas. 


r(s*/2+t)+ j ) 
:7*/2*s*/2+k) 


The  error  estimate  (38)  is  less  precise  than  (39);  however  it 
should  prove  easier  to  compute.  Also,  it  can  be  calculated  for  s  *  1. 
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To  use  (36)  to  approximate  G(r,  s,  a,  b;  k,  6,  x)  for  r  ■  1,  3,  5, 
s  ■  1,  3,  5,  by  replacing  the  infinite  series  with  a  partial  sum, 
we  must  of  course  be  able  to  compute  the  terms  of  the  series.  The 
evaluation  of  #[j+l,  r/2+s/2+k+j+l;  x(l-2a0)/(2a)]  can  be  accomplished 
as  described  above.  A  method  for  approximating  U[s/2+j,  r/2+s/2+j;  x(a*b)/(2ab) ] 
to  any  desired  accuracy  can  be  found  in  Section  3  under  Case  (iii) . 
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